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Aanleiding  tot  het  onderzoek 

In  Werkgebied  I  van  het  TNO  DO-onderzoeksprogramma  worden  rekenmodellen 
ontwikkeld  voor  constructiegeluidoverdracht  van  machines  naar  onderwater.  Een  van  de 
aspecten  die  daarbij  onderzocht  worden  is  het  effect  van  de  afmetingen  van  de  contact- 
vlakken  tussen  machine,  fundatie  en  scheepsbodem  op  de  geluidoverdracht.  In  een 
theoretische  studie  naar  lijnaanstoting  op  balk-achtige  constructies,  beschreven  in  TPD- 
HAG-RPT-930216,  is  aangetoond  dat  de  geluidoverdracht  ten  gevolge  van  stripaanstoting 
duidelijk  afwijkt  van  die  ten  gevolge  van  puntaanstoting,  waimeer  de  lengte  van  de  strip 
groter  of  gelijk  is  aan  de  helft  van  de  golflengte  in  de  balken.  Vanwege  het  praktische 
belang  voor  scheepsconstructies  is  het  onderzoek  voortgezet  met  een  studie  naar  de 
aanstoting  van  verstijfde  plaat-constructies  via  grote  contactvlakken. 

Doel  van  het  onderzoek 

Het  modelleren  van  aanstoting  van  een  plaat  met  een  lijnvormige  verstijver  via  grote 
contactvlakken.  Experimenteel  valideren  van  de  gevonden  modellen. 

Werkplan 

De  verdeelde  aanstoting  zal  worden  gemodelleerd  door  het  contactvlak  te  verdelen  in  een 
discreet  aantal  punten.  De  totale  energieoverdracht  wordt  vervolgens  bepaald  door  de 
ingangs-  en  overdrachtsadmittanties  en  de  krachtverdeling  voor  deze  punten  in  rekening 
te  brengen.  De  resultaten  van  modelberekeningen  zullen  worden  gevalideerd  aan  de  hand 
van  laboratoriumexperimenten. 

Conclusies 

Het  blijkt  gunsdg  te  zijn  de  aanstootkracht  onder  een  machine  te  verdelen  over  een 
contactvlak  met  tenminste  66n  dimensie  groter  dan  de  golflengte  in  de  scheepsconstructie. 
In  dit  rapport  worden  theoretische  en  experimentele  methoden  beschreven  waarmee  het 
effect  van  de  verdeelde  aanstoting  kan  worden  gekwantificeerd.  Twee  verschillende 
meetmethoden,  directe  admittantiemetingen  en  constructiegeluidintensiteitmetingen  op  de 
ontvangende  plaat,  leiden  tot  dezelfde  resultaten. 
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1  Introduction 


The  work  presented  in  this  report  has  been  performed  within  the  framework  of  the  TNO 
Defence  Research  programme  in  the  field  of  structural  acoustics.  This  programme  deals 
with  the  development  of  calculation  tools  for  the  prediction  of  the  effect  of  the  design  of 
the  seating  structure  on  the  vibratory  power  supplied  to  ship  structures  by  machinery. 

One  of  the  aspects  that  govern  the  structure-borne  noise  transmission  to  the  ship  structure 
is  the  size  of  the  connections  between  the  machine  or  its  seating  and  the  ship  structure. 
Conventionally,  noise  transmission  problems  are  treated  under  the  assumption  of  point¬ 
like  connections  between  source  and  receiver,  but  this  assumption  is  frequently  violated 
in  engineering  practice.  Thus  it  is  important  to  examine  the  effects  of  distributed 
excitation  via  large  contact  areas. 

An  introductory  study  of  surface  excitation  on  beam-  and  frame-like  systems  has  been 
presented  in  [1].  This  study  of  strip-excited  beams  and  frames  clearly  demonstrates  that 
effects  of  distributing  the  excitation  over  a  large  contact  area  can  be  both  beneficial  and 
detrimental  for  the  reduction  of  active  power  transmission  at  frequencies  where  the  strip 
length  is  of  the  order  of  magnitude  of  or  larger  than  the  governing  wavelength  in  the 
receiving  structure.  Similar  results  have  been  found  for  strip-excitation  on  an  infinitely 
large  homogeneous  plate  [2].  However,  the  effect  of  strip  excitation  in  the  vicinity  of  a 
structural  discontinuity,  like  a  stiffener,  on  a  plate-like  structure  remains  to  be 
investigated.  Therefore,  this  reports  describes  the  results  of  a  combined  theoretical  and 
experimental  study  of  strip  excitation  on  an  infinitely  large  homogeneous  plate  with  a 
single  line-stiffener. 


2 


Theoretical  background 


This  chapter  reviews  the  basic  definitions  and  describes  the  approach  chosen  for  the 
theoretical  study. 

2-1  Surface  mobility 

It  is  convenient  to  handle  structural  acoustic  power  transmission  by  means  of  mobility 
theory.  The  concept  of  surface  mobility  is  introduced  to  treat  the  variation  of  the  stress 
and  velocity  distributions  over  a  large  contact  area.  Far  a  flat  contact  area  S,  of  arbitrary 
shape,  in  the  xy-plane,  the  stress  distribution  c(x,y)  and  the  corresponding  velocity 
distribution  v(x,y)  yield  the  transmitted  complex  power  P  : 

P  =  L  Jo*(A:,y)v(x,y)dxdy  (1) 

The  total  force  F  is  obtained  by  integration  of  the  stress  distribution  over  the  contact 
area: 
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F  =  jCT(;c,)')dxd)^ 
s 


(2) 


The  velocity  distribution  is  related  to  the  stress  distribution  via  the  transfer  mobility 
distributions  '^(x,y\xf ,y' )  of  source  (subscript  's’)  and  receiver  (subscript  V): 

v(;t,y)  =  ^[y,(.x,y\x',y')  +  y  ^{x,y\x'  ,y')^G{x'  ,y')6x'  6iy'  (3) 

5 


In  analogy  with  the  definition  of  point  mobility,  the  surface  mobility  may  be  defined 
on  a  power  basis: 


y5  = 


P 


(4) 


or: 


s  _  s 


= 


jCT*(A:,y)|J[Y/A:,yljc',yO  +  y  ^{x,y\x' ,y')^<s{x' ,y')Ax' Ay'XdxAy 


Ja(j:,y)dA:dy 


(5) 


Hence  the  surface  mobility  is  weighted  by  the  stress  distribution  and  by  the  mobility 
distributions  of  source  and  receiver  in  the  contact  area. 


2.2 


Strip  mobility 


As  argued  in  [1]  it  is  worthwhile  to  study  strip  excitation,  i.e.  one  dimension  of  the 
contact  area  is  smaller  than  the  wavelength  while  the  other  can  be  large,  foremost 
because  of  an  enhanced  physical  lucidity.  The  problem  can  be  further  simplified  by 
describing  the  excitation  in  terms  of  a  stress  distribution  with  a  negligibly  small  source 
mobility  distribution  y,.  Any  source  can  thus  be  described,  by  selecting  the  right  stress 
distribution,  taking  into  accoxmt  the  actual  receiver  mobility.  Then  eq.(5),  for  a  strip  of 
length  L  along  the  x-axis,  simplifies  to: 


Y  — 


in.  in 

j'0*(x)|  J  y ^{,x\x')o{x')dx'\<i.x 

-Trt  ^ 


in 

2 

j  a(x)dx 

-in 


(6) 
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2.3  Effective  mobility 


The  effective  mobility  distribution  7*^  is  defined  as  the  resulting  mobility  distribution  on 
the  strip,  under  influence  of  the  complete  strip  excitation: 


y‘^{x)  = 


in. 

J  y ,(,x\x')<5(^x')(5.x' 

-in 


a(x) 


(7) 


Hence,  the  strip  mobility  is  related  to  the  effective  mobility  via: 


Y  strip 


La 

jy"^(x)  |cr(x)  pdx 

-in 

In  2 

J  CT(x)dX 
-in 


(8) 


2.4  Discretization 


For  the  beam-like  receiving  structure  there  exist  closed  form  solutions  of  the  above 
integral  equations  (see  [1]).  This  is  not  the  case  for  plate-like  receiving  structures  and 
certainly  not  for  a  plate  with  a  line  stiffener.  Therefore  the  integral  equations  for  the 
effective  mobility,  eq.(7),  and  for  the  strip  mobihty,  eq.(8),  are  solved  numerically  by 
means  of  discretization  of  the  stress  and  mobility  distributions.  The  strip  is  divided  in  N 
sections  of  equal  length  <ilj=UN.  In  each  section  n  the  stress  and  the  mobihty  are 
assumed  to  be  distributed  uniformly,  resulting  in  discrete  forces  F„  and  transfer  mobilities 
Ymn-  Hence  the  expression  for  the  effective  mobility  distribution,  eq.(7),  reduces  to: 


y'it  _  m-1 


Y  y  F 

•  nm  m 


(9) 


where  F„=C!(xJ6L.  The  strip  mobility,  eq.(8),  becomes: 


Y  strip 


N 


FP 


(10) 


where  F=IF^ ,  see  eq.(2). 

The  optimum  discretization  for  the  transfer  mobihty  distribution  is  governed  by  the 
wavelength  in  the  receiving  stracture.  Six  sections  per  structural  wavelength  should  be 
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sufficient.  Additionally,  the  disaetization  depends  upon  the  shape  of  the  stress 
distribution.  Smaller  sections  may  be  required  if  the  stress  distribution  exhibits  sharp 
peaks. 

2.5  Transfer  mobilities  on  a  line-stiffened  plate 

The  expressions  for  the  transfer  mobilities  on  an  infinitely  large,  homogeneous,  thin  plate 
with  or  without  a  single  line  stiffener  are  derived  in  [3]  and  experimentally  validated  in 
[4].  The  stiffener  is  modelled  by  means  of  a  rigid  line  support  along  the  y-axis,  i.e.  the 
plate  mobility  on  the  stiffener  is  zero.  The  expressions  are  given  in  Appendix  A. 


3  Numerical  modelling 


3,1  Strip  excitation  on  a  homogeneous  plate 

The  effect  of  strip  excitation  on  an  infinitely  large  homogeneous  thin  plate  in  flexure 
serves  as  a  reference,  to  be  able  to  separate  the  effects  of  distributed  excitation  from  the 
additional  effects  due  to  the  introduction  of  a  line  stiffener.  The  problem  is  studied  in  the 
range  of  Helmholtz-numbers  kj.,  based  on  the  flexural  wavenumber  and  strip  length  L, 
between  0.1  and  10.  Two  idealized  excitation  conditions,  a  uniform  force  distribution  and 
a  uniform  velocity  distribution  prescribed  along  the  strip,  are  examined. 

Uniform  force  distribution  (soft  strip) 

The  highest  Helmholtz-number  considered  (^^^10)  corresponds  with  the  shortest 
wavelength  (X,^=10/27rL).  Applying  the  criteriura  of  six  elements  per  wavelength  (see 
section  2.4)  leads  to  a  discretization  of  the  strip  into  A^=10  elements.  Figure  1  shows  the 
normalized  effective  mobility  {yfiNT)  under  the  strip,  for  Helmholtz-numbers  kbL=l, 

2,  5  and  10.  T”  is  the  input  mobility  for  point  excitation  on  an  infinite  plate,  x  is  the 
coordinate  under  the  strip  and  is  the  position  of  the  center  of  the  strip.  Because  of  the 
uniform  force  distribution,  this  mobility  corresponds  with  the  velocity  distribution  under 
the  strip.  It  is  seen  that,  with  increasing  Helmholtz-number,  the  real  part  of  the 
distribution  deviates  from  the  uniform  distribution  and  decreases.  The  resulting 
normalized  strip  mobility  (F'^Vr)  is  shown  in  Figure  2.  Above  Helmholtz-number  1 
the  real  part  of  the  strip  mobility  decreases,  which  indicates  that  less  power  is  transmitted 
via  the  strip  than  would  be  transmitted  via  point  excitation,  for  the  same  total  force. 

Uniform  velocity  distribution  (rigid  strip) 

The  discretization  applied  for  the  uniform  force  distribution  does  not  necessarily  apply 
for  the  uniform  velocity  distribution,  if  the  corresponding  stress  distribution  exhibits 
sharp  peaks.  The  stress  distribution  that  leads  to  a  uniform  velocity  distribution  follows 
from  inversion  of  the  mobility  matrix  for  the  discretized  strip,  eq.(9).  Figure  3 
shows  the  normalized  effective  impedance  distribution  (^„^  F”)  under  the  strip,  when  it 
is  discretized  into  A^=10  elements,  for  Helmholtz-numbers  2,  5  and  10.  Since  the 
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velocity  is  uniformly  distributed,  this  normalized  impedance  distribution  is  equal  to  the 
normalized  force  distribution.  The  force  appears  to  be  concentrated  at  the  ends  of  the 
strip.  Hence,  the  discretization  seems  to  be  insufficient.  Therefore  the  procedure  has  been 
repeated  for  N=20  and  A^=30  elements  respectively.  The  corresponding  force  distributions 
are  shown  in  Figure  4  and  Figure  5.  With  increasing  number  of  elements  the  force 
maxima  shift  towards  the  ends  of  the  strip.  Figure  6  shows  the  resulting  normalized  strip 
mobility  for  the  uniform  velocity  distribution,  with  iV=10,  20  and  30, 

respectively.  This  strip  mobility  is  calculated  on  the  basis  of  eq.(4)  and  the  calculated 
force  distribution.  The  strip  mobility  for  the  three  different  discretizations  do  not  differ 
significantly.  For  comparison,  the  strip  mobility  for  the  uniform  force  distribution 
(Figure  2)  is  added  in  the  same  figure.  Although  the  force  distributions  for  these  two 
cases  are  very  much  different,  the  strip  mobilities  show  a  comparable  behaviour.  The  real 
part  of  the  strip  mobility  for  the  rigid  strip  is  somewhat  lower  than  that  for  the  soft  strip, 
for  Helmholtz-numbers  between  circa  0.3  and  6. 

These  results  for  strip  excitation  on  a  homogeneous  plate  agree  with  those  reported  in 

[2]. 


3.2  Strip  excitation  on  a  line-stiffened  plate 

The  effect  of  introducing  a  line-stiffener  near  the  position  where  the  plate  is  excited  via  a 
strip  is  studied.  This  introduces  as  an  extra  parameter  the  position  of  the  center  of  the 
strip,  that  is  oriented  along  the  x-axis,  relative  to  the  fine-stiffener  along  the  y-axis  at  x=0. 
As  noted  in  section  2.5,  the  fine-stiffener  is  modelled  as  a  rigid  fine  support.  This  choice 
excludes  the  possibility  to  investigate  the  effects  of  excitation  by  a  rigid  strip,  that 
imposes  a  uniform  velocity  distribution,  when  the  strip  crosses  the  fine-stiffener,  i.e. 

Ix/U  <  1.  The  investigation  of  such  effects  would  require  an  adaptation  of  the  analytical 
model  for  the  mobilities  on  the  line-stiffened  plate  with  a  description  of  the  finite 
mobility  of  the  fine  stiffener. 

Figure  7  shows  the  velocity  distribution  under  the  strip  for  a  uniform  force  distribution, 
when  the  strip  is  centered  across  the  stiffener  (x/L=0).  The  velocity  at  the  stiffener 
equals  zero.  The  curves  in  Figure  7  do  not  show  this  because  they  are  drawn  between  the 
calculated  results  at  the  center  of  the  10  elements  into  which  the  strip  is  discretized.  At 
low  Helmholtz-numbers  the  effective  mobility  is  low,  since  the  stiffener  blocks  the  strip 
motion.  At  higher  Helmholtz-numbers,  i.e.  shorter  wavelengths,  where  the  plate  can 
deform  under  the  strip,  the  effective  mobility  increases.  The  same  behaviour  is  seen  in 
the  strip  mobility  (Figure  8).  Figure  9  shows  the  velocity  distribution  for  a  uniform  force 
distribution,  when  the  strip  just  touches  the  stiffener  (xJl=0.5).  In  this  case  the  strip 
starts  to  rotate  at  low  Helmholtz-numbers,  which  can  be  most  clearly  seen  for  The 
real  part  of  the  strip  mobility  (Figure  10)  exhibits  a  maximum  at  kjj=3.5.  Figure  11  and 
Figure  12  show  the  effective  mobility  and  the  strip  mobility  for  a  strip  that  does  not 
touch  the  stiffener  (xyl^l.O).  In  Figure  12  the  strip  mobilities  are  compared  for  a 
uniform  force  distribution  and  a  uniform  velocity  distribution.  In  the  calculation  for  the 
uniform  velocity  distribution  the  strip  is  discretized  in  30  elements,  to  capture  the  stress 
distribution  that  is  predominantly  concentrated  at  the  strip  ends,  as  has  been  discussed  in 
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section  3.1.  The  small  discontinuity  in  the  curves  for  the  uniform  velocity  distribution  at 
A:j2>1.5  is  caused  by  the  use  of  different  expressions  for  the  transfer  mobilities  for  low 
and  high  Helmholtz-numbers,  see  Appendix  A.  Small  deviations  in  the  transfer  mobilities 
are  amplified  in  the  matrix  inversion  that  is  applied  to  calculate  the  strip  mobility.  This 
effect  is  of  minor  importance  and  does  not  influence  the  general  conclusions.  Figure  13 
and  Figure  14  show  the  mobilities  for  a  strip  that  is  somewhat  further  removed  from  the 
stiffener  (x/[j=\.5).  It  can  be  seen  that  the  maxima  in  the  real  and  imaginary  parts  of  the 
strip  mobility  shift  towards  lower  Helmholtz-numbers  with  an  increasing  distance 
between  the  center  of  the  strip  and  the  stiffener. 

Figure  15  summarizes  the  effect  of  the  strip  and  the  line-stiffener  on  the  real  part  of  the 
strip  mobility,  i.e.  on  the  power  transmission.  In  this  case  the  strip  mobility  is  shown  for 
Helmholtz-numbers  up  to  100.  The  difference  between  the  strip  mobilities  for  a  uniform 
force  distribution  and  a  uniform  velocity  distribution  appears  to  be  unimportant.  As  a 
general  trend,  the  distributed  excitation  under  a  strip  results  in  a  decrease  with  a  factor  of 
2  (3  dB)  for  each  doubling  of  the  Helmholtz-number  k^L,  starting  at  /fcfrL=2.  The  presence 
of  the  line-stiffener  causes  a  sharp  decrease  of  the  strip  mobility  towards  low  Helmholtz- 
numbers.  The  value  of  kj^  below  which  the  mobility  starts  to  drop  off  decreases  with 
increasing  distance  between  the  center  of  the  strip  and  the  stiffener.  For  a  uniform  force 
or  velocity  distribution  under  the  strip,  the  overall  strip  mobility  is  the  lowest  when  the 
strip  is  centered  across  the  stiffener. 

The  results  of  these  numerical  simulations  agree  very  well  with  what  has  been  found  for 
distributed  excitation  on  beam-  and  fi'ame-like  systems  [1]. 


3.3  Area  excitation 

In  addition  to  the  strip  mobilities  discussed,  the  same  procedure  can  be  applied  to 
calculate  the  effect  of  distributed  excitation  of  a  two-dimensional  area  on  a  plate.  Some 
preliminary  calculations  have  been  performed  to  investigate  the  effect  of  extending  the 
width  of  the  strip.  A  uniform  force  distribution  over  a  square  area  of  dimensions  LxL  is 
considered,  both  on  a  homogeneous  plate  and  on  a  line-stiffened  plate  (for  x=0).  The 
area  is  discretized  in  1(30  equal  square  elements.  Figure  16  shows  the  resulting  real  parts 
of  the  surface  mobility  for  Helmholtz-numbers  k^^  up  to  50.  The  real  parts  of  the 
corresponding  strip  mobilities  are  displayed  as  a  reference  in  the  same  figure.  It  can  be 
seen  that  the  mobilities  for  the  strip  and  the  square  area  are  almost  equal  at  low 
Helmholtz-numbers,  but  that  at  higher  Helmholtz-numbers  the  surface  mobilities  for  the 
square  area  excitation  decrease  much  faster  with  increasing  Helmholtz-number  than  the 
strip  mobilities.  One  may  tentatively  conclude  that  it  is  advantageous  to  distribute  the 
total  force  of  a  machine  over  a  larger  contact  area,  because  this  reduces  the  energy 
Iransnussion  at  higher  frequencies,  although  the  present  calculations  only  support  this  for 
a  uniform  conphase  force  distribution.  The  investigation  of  the  effects  of  two-dimensional 
area  excitation  for  different  excitation  conditions  falls  without  the  scope  of  the  present 
smdy. 
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4  Experiments 


The  results  of  the  numerical  simulations  of  the  strip  mobility  that  have  been  described  in 
sections  3.1  and  3.2  are  validated  by  means  of  laboratory  experiments. 

4.1  Experimental  arrangement 

The  experiments  have  been  performed  on  the  same  line  stiffened  perspex  plate  that  has 
been  used  in  [4],  see  Figure  17.  It  is  a  20  mm  thick  perspex  plate  of  2  x  3  m,  hanging 
vertically,  on  rubber  strings  in  order  to  decouple  it  from  vibrations  of  the  building.  The 
line-stiffener  is  a  perspex  beam  of  thickness  10  mm  and  height  150  mm.  The  plate  is 
excited  by  a  shaker,  via  an  aluminum  strip  of  112  mm  long,  10  nun  wide  and  40  mm 
high.  In  one  of  the  experiments  it  has  been  attempted  to  create  a  uniform  force 
distribution  under  the  strip,  by  applying  a  thin  rubber  layer  between  the  strip  and  the 
plate,  see  Figure  18.  If  the  transfer  impedance  of  the  rubber  layer  is  much  smaller  than 
the  input  impedances  of  the  strip  and  the  plate,  then  it  converts  the  uniform  velocity  of 
the  rigid  strip  into  a  uniform  force  distribution  on  the  plate.  The  driving  force  on  top  of 
the  aluminum  strip  is  measured  by  a  force  transducer,  the  acceleration  of  the  strip  by  a 
set  of  two  small  accelerometers,  the  signals  of  which  are  added  to  suppress  the  rotational 
motion  of  the  strip.  This  rotational  motion  is  expected  to  be  negligibly  small,  because  the 
shaker  virtually  prevents  rotation. 

4.2  Characterization  of  the  rubber  layer 

As  argued  in  section  4.1,  the  force  distribution  on  the  plate  is  governed  by  the  velocity  of 
the  aluminum  strip  and  the  transfer  impedance  of  the  rubber  layer.  The  latter  has  been 
determined  experimentally,  via  the  method  that  is  used  at  TNO-TPD  to  characterize 
flexible  mountings  [5].  The  rubber  layer  is  mounted  on  a  cylindrical  steel  block  of  mass 
Af=18  kg.  The  aluminum  strip  is  used  as  top  mass.  It  is  driven  by  a  shaker,  like  in  the 
arrangement  shown  in  Figure  18.  The  layer  stiffness  K  and  loss  factor  Tj  are  calculated 
from  the  measured  complex  transfer  function  between  the  accelerations  above  (a,)  and 
below  (flj,)  the  layer: 


A:(1  +IT1)  =  -Moo^ _ (11) 

1  -  («./«,) 

The  results  are  shown  in  Figure  19.  The  approximate  stiffness  of  the  rubber  layer  is 
K=9  MN/m.  The  loss  factor  appears  to  be  frequency  dependent:  r|=0.007V/,  where /is 
the  frequency  in  Hz. The  deviations  from  pure  stiffness  behaviour  at  frequencies  above 
circa  800  Hz  indicate  the  presence  of  resonances  in  the  strip. 
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4.3  Measuring  the  strip  mobility 

It  is  not  straightforward  to  determine  the  strip  mobility  from  the  experiment.  The 
input  mobility  F""“  that  is  measured  on  top  of  the  aluminum  strip  has  to  be  corrected  for 
the  transmission  through  the  system  between  the  force  transducer  ad  the  plate  (i.e.  the 
connecting  rod,  the  aluminum  strip  and,  where  applicable,  the  rubber  layer).  This 
coirection  may  be  performed  by  modelling  the  aluminum  strip  with  the  connecting  rod  as 
a  mass  (m=0.21  kg,  determined  by  weighing)  and  the  rubber  layer  as  a  spring 
(A?=/!r(l+iTi),  see  section  4.2).  The  strip  mobility  is  then  calculated  from: 

ystrip  ^  -ioyfk  +  (1 

1  - 

An  alternative  way  to  determine  the  real  part  of  the  strip  mobility  is  to  measure  the 
power  transmitted  to  the  plate  and  divide  this  power  by  the  excitation  force  on  the  plate, 
i.e.  the  measured  excitation  force  corrected  for  the  acceleration  of  the  mass  m. 

4.4  Power  flow  measurements 

The  real  part  of  the  strip  mobility  can  be  determined  by  measuring  the  power  that  is 
transmitted  to  the  plate.  This  has  only  been  done  for  a  strip  position  that  is  remote  from 
the  stiffener.  In  that  case  the  power  transmission  to  the  plate  can  be  determined  from 
structural  intensity  measurements  on  the  plate,  around  the  strip.  The  validity  of  this 
method  has  been  checked  first,  for  point  excitation  in  a  position  at  a  distance  x,=400  mm 
from  the  stiffener,  close  to  the  center  of  the  plate  in  the  y-direction.  The  structural 
intensity  is  measured  in  Afc=10  positions  on  a  half  circle  of  radius  R=200  mm  around  the 
excitation  position,  see  Figure  20.  At  each  position  n  the  structural  intensity  I/n)  is 
determined  from  the  imaginary  part  of  the  cross-spectrum  of  two  accelerometer 
signals  [6]; 


IM) 


1  K 

8y"co^  sin(k^A) 


Im{G^^,(«)} 


(13) 


where  is  the  flexural  wavenumber,  Y“ =0.37"10'^  m/Ns  the  driving  point  mobility  of 
an  infinite  plate  (20  mm,  perspex)  and  A=32.5  mm  the  separation  of  the  accelerometers. 
The  radius  R  is  sleeted  to  maximize  the  distance  between  the  intensity  measurement 
positions  and  both  the  excitation  position  and  the  line  stiffener  (R=xJ2).  The  two-channel 
method  is  valid  for  flexural  wavelengths  smaller  than  2R,  i.e.  for  frequencies  above  circa 
500  Hz  (ki^~2n).  The  separation  of  the  accelerometers  is  chosen  such  that  it  is  smaller 
than  one  fourth  of  a  flexural  wavelength  at  the  highest  frequency  (5(X)0  Hz).  The  total 
power  flow  is  calculated  from  the  measured  intensities,  assuming  that  the  power  flow 
through  the  half  circle  that  is  not  measured  is  equal  to  the  measured  power  flow,  because 
of  symmetry  across  the  x-axis: 
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P 


n=l 


(14) 


The  input  mobility  is  related  to  the  power  flow  via  eq.(4).  The  driving  force  F  is 
measured  at  the  excitation  point.  Figure  21  shows  the  real  part  of  the  input  mobility, 
normalized  with  T”,  for  point  excitation  at  a  distance  of  400  mm  from  the  line-stiffener. 
The  solid  hne  gives  the  directly  measured  mobility,  using  a  force  transducer  and  a  set  of 
two  accelerometers  on  the  indenter  at  the  excitation  position,  the  dashed  line  is  derived 
from  the  intensity  measurements.  The  results  of  boths  measurements  agree  very  well.  The 
input  mobility  at  this  relatively  large  distance  from  the  stiffener  is  almost  equal  to  the 
infinite  plate  driving  point  mobility  F°.  The  measured  cuves  deviates  from  a  straight  line 
at  YIY°=\  because  of  the  modal  behaviour  of  the  non-infinite  experimenal  plate.  This 
comparison  shows  that  both  methods  can  be  used  to  determine  the  power  transmission  to 
the  plate. 

4.5  Strip  mobility  measurements 

The  strip  mobility  has  been  measured  at  position  x=400  mm  (L=112  mm,  i.e.  x/L=3.6), 
for  two  different  stress  distributions:  a  ’rigid’  strip,  i.e.  the  aluminum  block  mounted 
directly  on  the  plate,  and  a  ’soft’  strip,  with  a  rubber  layer  inserted  between  the 
aluminum  block  and  the  plate.  The  input  mobility  is  measured  on  top  of  the  block  and 
corrected  for  the  mass  and  stiffiiess  of  the  block  and  the  rubber  layer,  to  obtain  the  strip 
mobility.  The  real  part  of  the  strip  mobility  is  also  determined  from  intensity 
measurements.  Figure  22  shows  the  results  of  both  mehods  for  the  ’rigid’  strip.  Figure  23 
for  the  ’soft’  strip.  Again  (see  section  4.4)  the  agreement  between  the  results  of  both 
measurement  methods  is  excellent.  The  general  trend  of  the  results  is  seen  to  follow  the 
simple  model  that  has  been  found  from  the  numerical  modelling  (see  chapter  3),  for  strip 
excitation  remote  from  stiffeners: 


Re{T""'’}  ==  T” 

for 

k^L  <  2 

Re{r""M  = 

KL 

for 

k,L>2 

The  dotted  lines  in  Figine  22  and  Figure  23  give  the  mobility  as  measured  on  top  of  the 
aluminum  strip.  At  high  Helmholtz-numbers  this  mobility  is  lower  than  the  strip  mobility 
due  to  the  acceleration  of  the  mass  of  the  block.  Figure  23  shows  the  damped  resonance 
of  the  block  on  the  rubber  layer  at  about  1000  Hz  (kip=2.4).  Figure  24  shows  the  real 
and  imaginary  parts  of  the  strip  mobihty,  derived  from  the  input  mobility  measurement, 
together  with  the  results  of  theoretical  predictions,  according  to  the  numerical  modeling 
that  has  been  described  in  chapter  3.  The  average  trend  of  the  measured  imaginary  parts 
of  the  strip  mobility  starts  to  deviate  from  the  theoretical  predictions  above  kjL=3.  An 
explanation  for  this  deviation  has  not  been  found.  The  modal  behaviour  of  the  receiving 
plate  largely  masks  the  small  differences  that  might  be  expected  due  to  the  different 
stress  distributions,  although  one  might  tentatively  conclude  that  the  measured  mobilities 


Ongerubriceerd 


TNO- report 


Ongerubriceerd 


TPD-HAG-RPT-9601 32 


confirm  the  prediction  that  the  mobility  for  the  soft  strip  (’uniform  force’)  is  somewhat 
higher  than  that  for  the  rigid  strip  (’uniform  velocity’)  at  Helmholtz-numbers  around 
k^2. 

Figure  25  shows  the  measured  mobilities  for  positions  on  the  stiffener.  The  driving  point 
mobihty  on  the  stiffener  is  compared  with  strip  mobilities  for  a  rigid  strip  at  x/L=0  and 
x/L=021.  The  strip  mobilities  are  almost  equal  for  both  positions.  They  do  not  show  the 
sharp  decrease  of  the  mobility  towards  low  Helmholtz-numbers  that  is  found  in  the 
numerical  modelhng  (see  Figure  15).  The  reason  for  this  is  clear  from  the  comparison 
with  the  measured  point  mobility  on  the  stiffener.  In  the  numerical  model  the  line 
stiffener  is  modelled  as  a  rigid  support,  hence  this  point  mobility  is  infinitely  small.  In 
the  experiment  the  line  stiffener  has  a  finite  stiffness,  so  that  the  strip  mobility  at  low 
Helmholtz-numbers  is  governed  by  the  input  mobility  of  the  stiffener.  The  driving  point 
mobility  increases  towards  kjj=5  where  it  appproaches  a  resonance  for  longitudinal  (’in¬ 
plane’)  waves  in  the  height  of  the  stiffener.  This  behaviour  is  not  seen  in  the  strip 
mobilities,  because  the  strip  sees  the  mobility  of  the  stiffener  and  the  strip  mobility  on 
the  plate  in  parallel,  so  that  the  lowest  mobility  predominates.  At  k^->5  the  strip  mobility 
is  lower  than  the  point  mobility  on  the  beam,  as  may  been  seen  from  comparison  with 
Figure  24.  As  stated  in  the  beginning  of  section  3.2,  the  present  model  for  strip  mobility 
on  a  line-stiffened  plate  would  require  an  adaptation  to  include  the  effect  of  the  mobility 
of  the  stiffener. 
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Conclusions  and  recommendations 


The  distribution  of  the  excitation  force  from  a  machine  on  the  ship  structure  over  a  large 
contact  area  is  beneficial  for  the  reduction  of  noise  transmission  at  higher  frequencies. 
The  criterium  for  the  ocurrence  of  this  effect  is  that  the  contact  area  should  have  at  least 
one  dimension  larger  than  the  governing  wavelength  in  the  ship  structure. 

This  general  conclusion  may  be  drawn  from  the  results  of  the  theoretical  and 
experimental  investigations  described  in  this  report.  The  effect  can  be  quantified  in  terms 
of  a  ’surface  mobility’,  that  is  defined  as  the  ratio  of  the  total  transmitted  power  via  the 
contact  area  and  the  total  excitation  force  in  the  contact  area.  Numerical  predictions  can 
be  made  by  applying  previously  derived  analytical  expressions  for  driving  point  and 
transfer  mobilities  on  plate  fields  with  and  without  line-stiffener  in  combination  with  a 
discription  of  the  force  distribution  in  the  contact  area.  From  numerical  simulations  it 
follows  that  the  differences  in  the  power  transmission  are  small  for  two  idealized 
excitation  conditions,  a  uniform  conphase  force  distribution  and  a  uniform  conphase 
velocity  distribution,  in  spite  of  large  differences  in  the  corresponding  spatial  force 
distributions.  This  conclusion  is  validated  in  a  series  of  laboratory  experiments. 

Two  different  experimental  techniques  have  been  applied  to  determine  the  strip  mobility 
for  excitation  of  a  line-stiffened  plate  by  a  rigid  block,  that  is  either  mounted  directly  on 
the  plate,  or  mounted  resiliently  via  a  thin  rubber  layer.  Direct  input  mobility 


Ongerubriceerd 


TNOreport 


Ongerubriceerd 


TPD-HAG-RPT-960132 


measurements  on  top  of  the  block,  corrected  for  the  transmission  via  the  block  and  its 
mounting,  and  structural  intensity  measurements  on  the  receiving  plate  yield  the  same 
results  for  the  strip  mobility. 

An  important  conclusion  of  this  study  is  the  fact  that  the  concept  of  a  ’surface  mobility’ 
is  valid.  This  offers  the  possibility  to  implement  the  structure-borne  sound  transmission 
via  large  contact  areas  in  a  prediction  method  in  terms  of  a  single  mobility.  The 
prediction  of  the  effects  of  surface  excitation  across  a  stiffener  on  a  plate  would  require 
an  extension  of  the  available  model  for  the  rigidly  supported  plate  with  a  description  of 
the  finite  mobility  of  the  line-stiffener.  Furthermore  it  is  highly  relevant  to  acquire 
information  about  the  actual  force  distributions  under  the  mounts  of  shipboard  machinery. 


Delft,  23  September  1996  TNO  Institute  of  Applied  Physics 


dr.ir.  C.A.F.  de  Jong,  ir.  P.P.  Kooyman,  ing.  F.G.P.  van  der  Knaap 
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strip  on  homogeneous  plate:  unifonn  force  distribution 


Figure  1:  real  (top)  and  imaginary  (bottom)  part  of  the  normalized  effective  mobility 
distribution  under  a  strip  with  uniform  force  distribution  on  an  infinite  homogeneous 
plate. 


strip  on  homogeneous  plate:  uniform  force  distribution 


Figure  2:  real  (top)  and  imaginary  (bottom)  part  of  the  normalized  strip  mobility  for  a 
strip  with  uniform  force  distribution  on  an  infinite  homogeneous  plate. 
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(x~xc)/L 


Figure  3:  real  (top)  and  imaginary  (bottom)  part  of  the  normalized  effective  impedance 
distributwn  under  a  strip  with  uniform  velocity  distribution  on  an  infinite  plate;  N=10. 


Figure  4:  real  (top)  and  imaginary  (bottom)  part  of  the  normalized  effective  impedance 
distribution  under  a  strip  with  uniform  velocity  distribution  on  an  infinite  plate;  N=20. 
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(x-xc)/L 


Figure  5:  real  (top)  and  imaginary  (bottom)  part  of  the  normalized  effective  impedance 
distribution  under  a  strip  with  uniform  velocity  distribution  on  an  infinite  plate;  N=30, 


kbL 


Figure  6:  real  (top)  and  imaginary  (bottom)  part  of  the  normalized  strip  mobility  for  a 
strip  with  uniform  velocity  on  an  infinite  plate  (N=10,  20  &  30)  and  for  a  uniform  force 
distribution. 
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Strip  on  line-stiffened  plate:  uniform  force  distribution;  xc/L=0 
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Figure  7:  real  (top)  and  imaginary  (bottom)  part  of  the  normalized  effective  mobility 
distribution  under  a  strip  with  uniform  force  distribution  on  a  line-stiffened  plate;  x/L=0, 


strip  on  line-stiffened  plate:  unifcon  force  distribution;  xc/L=0 


Figure  8:  real  (top)  and  imaginary  (bottom)  part  of  the  normalized  strip  mobility  for  a 
strip  with  uniform  force  distribution  on  a  line-stiffened  plate;  x/L=0. 
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strip  on  line-stiffened  plate:  uniform  force  distribution;  xc/L=0.5 
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Figure  9:  real  (top)  and  imaginary  (bottom)  part  of  the  normalized  effective  mobility 
distribution  under  a  strip  with  uniform  force  distribution  on  a  line-stiffened  plate; 
x/L=0.5. 


Strip  on  line-stiffened  plate:  uniform  force  distribution;  xc/L=0.5 


Figure  10:  real  (top)  and  imaginary  (bottom)  part  of  the  normalized  strip  mobility  for  a 
strip  with  uniform  force  distribution  on  a  line-stiffened  plate;  x/L=0.5. 
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(x-xc)/L 

Figure  11:  real  (top)  and  imaginary  (bottom)  part  of  the  normalized  effective  mobility 
distribution  under  a  strip  with  uniform  force  distribution  on  a  line-stiffened  plate; 
x/L=1.0. 


Figure  12:  real  (top)  and  imaginary  (bottom)  part  of  the  normalized  strip  mobility  for  , 
strip  with  uniform  velocity  or  force  distribution  on  a  line-stiffened  plate;  x/L=1.0. 
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Figure  13:  real  (top)  and  imaginary  (bottom)  part  of  the  normalized  effective  mobility 
distribution  under  a  strip  with  uniform  force  distribution  on  a  line-stiffened  plate; 
x/L=1.5. 


Strip  on  line-stiffened  plate:  uniform  velocity  distribution  (xc/L=1.5) 


Figure  14:  real  (top)  and  imaginary  (bottom)  part  of  the  normalized  strip  mobility  for  a 
strip  with  uniform  velocity  or  force  distribution  on  a  line-stiffened  plate;  x/L=L5, 
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strip  mobility 


Figure  15:  the  real  part  of  the  normalized  strip  mobility  on  a  homogeneous  and  a  line 
stiffened  plate,  in  dB,  compared  for  different  stress  distributions  under  the  strip. 


Figure  16:  the  real  part  of  the  normalized  mobility  for  uniform  force  distribution  on  a 
homogeneous  and  a  line-stiffened  plate,  in  dB;  for  excitation  via  a  strip  of  length  L  and 
via  a  square  area  (LxL) 
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Figure  1 7:  Sketch  of  the  perspex  plate  on  which  the  experiments  have  been  performed. 


Figure  18:  Photograph  of  the  aluminum  strip,  mounted  on  the  plate  via  a  rubber  layer. 

The  strip  is  driven  by  a  shaker,  via  a  force  transducer.  The  acceleration  of  the  strip  is 
measured  by  a  set  of  two  small  accelerometers. 
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rubber  layer 


100  200  500  Ik  2k  3k  4fc 
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Figure  19:  stiffness  (top)  and  loss  factor  (bottom)  of  the  rubber  layer  that  is  applied  for 
the  'uniform  force'  strip  excitation.  The  dashed  curves  are  fitted  to  the  measured  data. 


Ongerubriceerd 


TNO-report 


Ongerubriceerd 


TPD-HAG-RPT-960132 


25 


Figure  20:  photograph  of  the  power  flow  measurement  for  point  excitation,  using  a 
pendulum.  Structural  intensity  measurements  are  performed,  using  two  accelerometers, 
on  10  positions  on  a  half  circle. 


point  measurement  on  line-stiffened  plate  (xs={).4m) 


Figure  21:  Normalized  real  part  of  the  input  mobility  for  point  excitation  at  x=0.4  m, 
measured  directly  and  determined  from  intensity  measurements. 
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Figure  22:  normalized  real  part  of  the  mobility  for  a  rigid  strip  (x/L=3,6f  measured 
directly,  with  and  without  correction  for  the  strip  mass,  and  determined  from  intensity 
measurements. 
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Figure  23:  normalized  real  part  of  the  mobility  for  a  soft  strip  (x/L=3.6),  measured 
directly,  with  and  without  correction  for  the  strip  mass  and  stiffness,  and  determined 
from  intensity  measurements. 
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strip  on  line-stiffened  plate;  xc/L=3.6 


Figure  24:  real  (top)  and  imaginary  (bottom)  part  of  the  normalized  strip  mobility  at 
x/L=3.6.  Measurements  for  rigid  and  soft  strips  compared  with  theoretical  predictions. 
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Figure  25:  real  (top)  and  imaginary  (bottom)  part  of  the  normalized  mobility  for  point 
excitation  in  x/L=0  and  excitation  via  a  ’rigid’  strip  in  x/L=0  and  x/L=0.27. 
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Appendix  A:  Transfer  mobilities  on  an  infinite  homogeneous  thin  plate 
with  a  single  line  discontinuity 

This  Appendix  reviews  the  expressions  for  the  transfer  mobility  (see  [3], [4]  and  [7]).  The 
plate  lies  in  the  xy-plane.  The  distance  between  the  excitation  ("source":  s)  and  receiver 
(r)  positions  is  defined  as: 

Z)„  = 


and  the  distance  from  the  receiver  position  to  the  perpendicular  projection  of  the  source 
position  on  the  discontinuity  as: 


D,  =  )lx^  + 


(A.2) 


The  bending  wave  propagation  on  an  infinite  plate  is  governed  by  the  propagation 
function: 


r(,(*,r)  =  -  Ho^\-ik,r)  (A.3) 

where  kg  denotes  the  bending  wavenumber,  r  the  distance  and  the  complex  Hankel 
function  of  the  second  kind  of  order  n.  The  ’first  order’  propagation  function  appears 
when  the  propagation  function  is  derived  with  respect  to  its  argument: 

n,(A:,r)  =  =  H,^\k,r}  +  iH,^\-ik,r)  (A.4) 

dk^r  " 


For  the  transfer  mobilities  on  an  infinite  plate  with  a  single  line  discontinuity  three 
different  approximate  expressions  are  used  for  different  conditions,  that  are  determined 
by  the  source  and  receiver  positions,  relative  to  the  line  discontinuity  along  the  y-axis  at 
x=0,  and  by  the  Helmholtz-number  (explained  below)  : 

-  for  \  /  x^  \  >  1  and  ’high’  Helmholtz  numbers: 

“for  \  Xr  /  x^  \  >  I  and  ’low’  Helmholtz  numbers: 

mw  (A.6) 

-for  \  x, !  x,  \  <  \  (reciprocal  expression)  : 
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)  (A.7) 

The  first  part  of  eq.(A.5),  i.e.  gives  the  transfer  mobility  on  a  homogeneous 

plate,  without  line  stiffener.  The  second  part  of  this  equation  expresses  the  effect  of 
reflection  at  the  line-stiffener. 

It  appears  to  be  necessary  to  apply  the  ’low’  Helmholtz  number  expression,  because  the 
’high’  Helmholtz  number  expression  underestimates  the  transfer  mobilities  for  positions 
on  the  same  side  of  the  discontinuity  and  overestimates  the  transfer  mobilities  across  the 
discontinuity.  This  leads  to  erroneous  results  in  the  integration  that  is  performed  to 
calculate  the  effective  and  the  strip  mobilities.  The  transition  point  between  ’low’  and 
’high’  Helmholtz  numbers  is  chosen  where  the  spectra  of  the  two  approximations  cross 
or  where  they  start  to  diverge.  The  Helmholtz  number  at  which  this  occurs  is  different 
for  the  real  and  the  imaginary  part  of  the  transfer  mobilities.  The  transition  point  has 
been  determined  visually  from  a  series  of  plots  of  transfer  mobilities  for  different  source 
and  receiver  piositions,  see  figures  Al,  A2,  A3  and  A4.  An  appropriate  choice  for  the 
transition  point  appears  to  depend  upon  the  Helmholtz  number  based  on  the  distance 
from  the  receiver  position  to  the  perpendicular  projection  of  the  source  position  on  the 
discontinuity  : 


transition  point  at  =  0.8  for  Ref  T  _  } 

"  ^  (A.8) 

transition  point  at  k^D^  =  0.1  for  Im{ 

Note  that  reciprocal  expression  is  taken  if  the  receiver  position  is  closer  to  the 
discontinuity  than  the  source  position.  In  that  case  the  transition  point  is  based  on  the 
after  interchanging  the  source  and  receiver  positions. 
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Figure  Al:  Transfer  mobility  for  x/x=l,  according  to  eq.(A.5)  and  eq.(A.6).  The  domed 
line  indicates  the  chosen  transition  point 


xrAcs  =  -1 


Figure  A2:  Transfer  mobility  for  x/x,--l,  according  to  eq.(A.5)  and  eq.(A.6).  The  domed 
line  indicates  the  chosen  transition  point 
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Figure  A3:  Transfer  mobility  for  x/x=10,  according  to  eqfAJ)  and  eq,(A.6).  The  domed 
line  indicates  the  chosen  transition  point. 


Figure  A4:  Transfer  mobility  for  x/x  =-10,  according  to  eq.(A.5)  and  eq.(A,6).  The 
domed  line  indicates  the  chosen  transition  point. 
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